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I. INTRODUCTION 



The present paper reviews properties of a three-dimensional weakly diluted quenched Ising model (random Ising 
model: RIM - see figure |^) in the vicinity of critical point. For weak dilution by a non-magnetic component, i.e. far 
from the percolation threshold, the magnetic second-order phase transition still maintains in the RIM although it 
^ ■ possesses novel features in comparison with the pure d = 3 Ising model. Static critical exponents of a RIM ^ have 
O ] been a subject of a detailed experimental P-pOf, numerical pl|-pO[ and theoretical E^-^ analysis almost for three 
, decades. Recently, new data were obtained both in experimental measurements ||l5|-|20|| and Monte Carlo simulations 

p3|-p0[. Theoretical breakthrough occurred along several months in 1999-2000 when the perturbation expansion 
T— I ; series for RIM were extended from the 4th through 5th to the sixth order |66 6^. Therefore it seems for 

K*" ' us worth collecting together the bulk of available results concerning critical behaviour of the RIM. 
00 ^ In the present discussion, we will be mainly interested in the RIM critical exponents. An asymptotic critical exponent 
' x of a physical observable 0(t) is defined |^ asymptotically close to the critical point Tc- x = limT-_>o '"^fj^ , where 
T — {T — Tc)/Tc is the reduced distance to the critical point. For instance, the magnetic susceptibility x diverges as: 

o ■ 

X-r±|r|-^ r^O (1) 

o . 

■ where 7 is the susceptibility critical exponent, while r+ and r_ are critical amplitudes above and below the critical 
^ J point respectively. The power law of type (|l|) holds exactly in the asymptotic regime r — > 0. In this regime the critical 
^ ■ exponents and ratios of the critical amplitudes take constant values. According to the universality hypothesis they 
are defined by global variables only. For a short-range interaction the global variables are the space dimension and 
tensor characteristics of an order parameter. In the non-asymptotic regime the approach to criticality is described by 
Q ' non-universal effective critical exponents^ which are introduced to describe the behaviour of a quantity in a certain 
i temperature interval [ |68|j69[| . The susceptibility effective critical exponent 7^^,^. by definition reads: 



l.Ar)^--^. (2) 



dlnx(T) 

^ In the asymptotic limit r ^ the effective and asymptotic exponents coincide. In the intermediate region the 
behaviour is characterized by the so-called Wegner expansion |]70| : 

X ^ r±|rr^ (1 + T,^±\Tr + r2^±|Tp'^'' + ...), (3) 

where ri_-|-,r2,± are non-universal amplitudes, v is the correlation length critical exponent, and lo is the correction- 
to-scaling exponent. 

It is the art of a physicist which allows to discriminate between regimes (0)~(^ either performing experimental 
measurements or doing theoretical calculations and numerical simulations. Below we will review available data on 
effective and asymptotic critical exponents of the RIM. The paper is arranged as follows. In the next section || we 
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will give general considerations about an influence of a weak quenched disorder on a second order phase transition, 
formulate the model as well as outline the main features of its critical behaviour. In Section III we review experiments 



on weakly diluted uniaxial magnets, whereas computer Monte Carlo simulations of a RIM are analyzed in Section 
IV. In sections and we outline the method of renormalization group as the most fruitful theoretical tool in the 
study of the critical properties of RIM. We consider different renormalization sche mes and pay special attention to 



analysis of divergent series obtained within the approach. Conclusions (Section |VI|) complete the article 



II. WEAKLY DILUTED QUENCHED ISING MODEL 

The central questions one has to answer studying the influence of weak disorder on a magnetic second-order 
phase transitions are: do the critical exponents of a homogeneous magnet change under dilution by a non-magnetic 
component? And if this is the case, are the new exponents universal? Regarding the first question it was argued |7l| ] 
that if the heat capacity critical exponent a of the pure (undiluted) system is positive, i.e. the heat capacity diverges 
at the critical point, then a quenched disorder causes changes in the critical exponents. This statement is known as 
Harris criterion. Later, for a large class of d-dimensional disordered systems it was proven that the correlation length 
critical exponent v must satisfy the bound v >2/d [Q. Both statements focus attention towards studies of the d — i 
Ising model, where the typical numerical values of the above exponents, together with the magnetic susceptibility and 
the order parameter critical exponents in the pure case read |73| : 

a = 0.109 ± 0.004 > 0, = 0.6304 ± 0.0013 < 2/3, 7 ^ 1.2396 ± 0.0013, (3 0.3258 ± 0.0014. (4) 

According the inequalities of the criteria mentioned above for the case of the diluted Ising model new exponents are 
expected. 

In order to obtain precise values of critical exponents it is standard now to rely on renormalization group methods. 
In particular, the theoretical estimates given above (jj) were obtained on the basis of a deep analogy between long- 
distance properties of the Ising model in the neighbourhood of a second-order phase transition point and a field theory 
with an effective Landau-Ginzburg- Wilson Hamiltonian of the form: 

Hi.i„g(^) = j d^R[\ [\V^\'+ml^^] + (5) 

where ttiq is a bare mass proportional to distance to a critical point, tp — ^{R), uq are (bare) scalar field and coupling. 

One of the ways to introduce quenched dilution to the effective Hamiltonian (||) is to add to uiq the random- 
temperature like variable ^ — ip{R): 

n^icp) = J d'R{l [IV^P + (ml + ij)^^] + (g) 

Taking that ip obeys Gaussian distribution: 

P(V) = ^^cxp(-^V4ii;2) (7) 



with w"^ being a half of dispersion, and introducing n replicas [ [75[ of a model (^) in order to perform averaging over 
quenched disorder U& one ends up Em with a familiar effective Hamiltonian: 



In the limit n ^ field theory (||) describes critical properties of a RIM. Here, the bare coupling uq is positive (being 
proportional to uq) whereas the bare coupling vq is proportional to minus variance of the random variable ip and thus 
is negative. The last term in (|8|) is present only for non-zero dilution: it is directly responsible for the fluctuations 
effective interaction due to the presence of impurities. 

Often (e.g. in Monte Carlo simulations) the microscopic Hamiltonian of RIM is written in the following form: 

R,R' 
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where R runs over sites of simple cubic lattice, J is a short-range translationally invariant interaction between pairs of 
5 = ±1 "classical" Ising spins, cr is an occupation number equal either to or to 1. It is considered that geometrically 
the vacancies cr = are distributed independently according to the law 

P{cr) = (1 - p)6{cii) + p5{l - cr) (10) 

and fixed in different sites of the lattice (see figure |^). In ( p^ ) p < 1 is the concentration of occupied sites. Calculating 
free energy of a model (^ and using replica trick ]7^ ] to perform the configurational averaging of the logarithm of 
configuration-dependent partition function it is straightforward to show that one ends up with the effective Haniilto- 
nian of (||) type in the n ^ replica limit with uo ^ p and vq ~ p{p — 1). 

From the viewpoint of dynamics one can point to two opposite types of a disorder |^ . If a characteristic time of 
impurities dynamics is comparable to relaxation times in the pure system, impurity variables are treated identically 
to the "pure" dynamical variables and are a part of the disordered system phase space. The corresponding annealed 
disorder is a subject of special investigations and reveals trivial results at criticality. The so called Fisher 
renormalization flT^ states that when a heat capacity critical exponent of an undiluted system apure is positive, then 
the critical exponents x of an annealed system are determined by those of the corresponding pure one (xpure) by a 
simple renormalization of the form: 

X —n. 

(11) 



1 — a 1 — a 

pure pure 

This explains why prevailing interest is attracted by the quenched disorder when impurities can be considered as 
fixed and thus one needs to perform configurational average over an ensemble of disordered systems with different 
realization of the disorder. 

The Hamiltonian (^ represents critical properties of the problem (^ for small randomness. Alternatively, the 
scale-invariant fractal ramified cluster at the percolation threshold is a starting point of the strong disorder approach. 
Here, the field theoretical description starts with the effective Hamiltonian of the Potts model. A unified theory of 
random systems critical behaviour which would give both regimes of strong disorder and weak disorder as its limiting 
cases is still absent. 

The translation invariant Hamiltonian (H) presumes perturbative account of thermal fluctuations around a spatially 
homogeneous unique ground state. While it appears to hold for a pure system, for a random system there is a 
macroscopic number of spatially inhomogeneous ground states in a disorder dominated region. These correspond 
to local minimum solutions of a saddle-point equation for the effective Hamiltonian (^) [ |78|]79| ]. Physically the last 
corresponds to the so-called Griffiths phase caused by the existence of ferromagnetically ordered "islands" in the 
temperature interval between the critical temperatures of pure and random systems. The description of the phase at 
the critical point is provided by a replica-symmetry breaking Hamiltonian | |78| leading to non-trivial results | pl| , p2| . 
However recently a refined analysis of the problem brought about a stability of the critical behaviour of the weakly 
disordered systems with respect to replica symmetry breaking effects pSll. The theoretical results reviewed in the 
present discussion will be based on the replica symmetrical Hamiltoniari(^ (see section [v[). 

The effective Hamiltonians (|^) and (j8|) possess different global variables: although the space dimension is the same 
(d = 3) the symmetry and order parameter components number differ. Thus one may expect that by application 
of renormalization group approach they will lead to different critical behaviour. One of the central notions in the 
renormalization group formalism applied to the critical phenomena is the notion of the fixed point of the renormal- 
ization group transformation. If a fixed point exists and it can be reached from the initial values of the couplings 
it corresponds to a critical point of the system. Applying the renormalization group transformation to the effective 
Hamiltonian (^) and starting from positive values of u one reaches the stable fixed point u* which corresponds to the 
critical point of a c? = 3 pure Ising model and leads in particular to the results (^. The fixed point structure for 
the effective Hamiltonian (H) is sketched in figure 0. The earliest qualitative results about the structure of RIM fixed 



points appeared in mid-seventies [^-45|. Later analysis supported this picture [ p8pC| , |84| ] : indeed the fixed point I 
of the d = 3 pure Ising model appears to be unstable and a new stable fixed point R exists (see figure |^) . Thus the 
general answer of the renormalization group analysis of the RIM supports non-perturbative results of Refs | |7l| , [72| : 
RIM critical behaviour is governed by critical exponents different from those of the pure Ising model. 

In the subsequent chapters we will have a look how this statement was made clear in experimental, Monte Carlo 
and theoretical studies. 

III. EXPERIMENTAL STUDY 

The Ising model is represented in experiments on magnets by crystalline difluoride of a transition metal, normally 
iron or manganese. In FeF2 which possesses a rutile crystalline structure with a = 4.697^1 and c = 3.309A at room 
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temperature the spins of metal ions align along c axis in such a way that spins of the body-centered ions are 
oriented in opposite directions to those of the corner ions. The system can be well described by a Heisenberg-like 
spin S = 2 Hamiltonian with quadratic anisotropic single ion terms and other interactions less than 6% of interspin 
exchange force constant J = 0.45 meV, as follows from the spin-wave dispersion relations on the basis of neutron 
scattering measurements [ |8^ . Dominant intersublattice exchange interaction and high anisotropy makes FeF2 a very 
good experimental realization of an Ising antiferromagnet, where the order parameter is the sublattice magnetization. 
Experimental studies confirm pure Ising model critical behaviour for the reduced temperature range |r| < 10""'^ 
(r = {T — Tis[)/T^) around Neel temperature Tn- Another example is MnF2 with similar crystalline structure but 
much weaker single ion anisotropy than FeF2, however the experimental study showed the substance to belong to the 
Ising universality class ||87|-^ as well. 

A material which corresponds to RIM can be obtained as a crystalline mixture of two compounds on the basis of 
"pure" Ising model matrix (see Table ^. A corresponding site-diluted uniaxial alloy can be prepared by substitution 
of FeF2 (MnF2) with non-magnetic isomorph ZnF2. Experiments on the critical behaviour of random systems are 
extremely sensitive to the sample quality. Asymptotic critical behaviour is observable only very close to Tn and 
macroscopic non-statistical gradients in concentration cause variation of Tn through the sample smearing out the 
sharp transition. Therefore, to provide a satisfactory good realization of random substitutional disorder of magnetic 
ions (Fe+^, Mn+^) by non-magnetic ones (Zn+^), mixed crystals FepZni_pF2 as well as MnpZni_pF2 should be grown 
with very high crystalline accuracy of the mosaic pattern, chemical homogeneity and especially with small impurity 
concentration gradients. The last can be achieved by choosing a concentration of impurities such that pl| dT^/dp = 
or by a geometry which takes into account that usually gradients are parallel to the growth axis. Early experimental 
studies on the disordered crystals critical behaviour |92 9^ proved the crucial role of a sample quality and high 
quality samples allowed to observe sharp phase transition and to measure the dependence of Neel temperature Tpf 
on p via linear birefringence method |94[ . This provided a basis for universal critical properties measurements. The 
earliest experimental study of the critical exponents governing the sharp transition at a weak quenched dilution was 
the nuclear magnetic resonance measurement of magnetization in Mno.864Fo.i36F2 The value of magnetization 
exponent f3 was found to differ strongly from that in non-diluted sample (see Table |l for details) . 

In a few years the study of Rcf. H was corroborated by nuclear scattering (NS) measurements in a two-axis 
spectrometer configuration of the staggered susceptibility and the correlation length in a crystal on the iron 
basis with concentration p — 0.5 in the high dilution regime. In the experiment the smearing effect was eliminated by 
masking the crystal to expose only a small homogeneous region, while the relative temperature control to about O.OIK 
together with the substance Neel temperature 42.50if gave the accuracy S\t\ ~ 5 x lO"'' for reduced temperature. 
The data for the inverse correlation length and the inverse staggered susceptibility were fit to power laws with critical 
amplitudes and exponents as free parameters. The power law was shown to be satisfied well within lO^'^ 2 x 10"^ 
with the same critical exponents (see Table ^ below and above Tc and differing strongly from those of the pure Ising 
model. To support the result by alternative experimental methods caloric properties were measured. Utilizing the 
fact that for transition-metal difluorides the temperature derivative of the linear birefringence is proportional to the 
magnetic specific heat within ||9^ , |95|] |r| < lO"'^, measurements of the critical exponent a were performed for a sample 
Feo.6Zno.4F2 with Tn = 47.05. In order to minimize the effect of concentration gradients the laser beam was oriented 
perpendicularly to the concentration gradient. The numerical value of the critical exponent was extracted by fitting 
the data to the temperature integral of specific heat scaling function taking into account correction-to-scaling terms 
(i.e. the first two terms in the expansion (||)). In contrast to the pure Ising case, the data obtained yielded a < (see 
Table |). Within no region in reduced temperature r any evidence of pure Ising behaviour or a crossover from pure 
to random fixed point critical behaviour was found. So one concluded that the crossover is either outside the critical 
region or it is too slow. The data obtained in Ref. on the basis of two different substances and quite different 
experimental procedures proofed that the scaling relation dv — 2 — a holds for them and thus supported strongly 
scaling in dilute systems. 

With a typical energy of the neutron beam about 10 meV, the NS measurements appeared to be one of the most 
useful method to study criticality, however they turned out to be in particular very sensitive to the sample quality since 
typical sample sizes are less than several millimeters. The NS measurements contain contributions from transverse 
and longitudinal spin fluctuations, which have to be separated in a detailed data analysis. This explains why in 
the neutron scattering measurements of the Ref. Q which were performed on FepZni_pF2 with concentration 0.46 
on the similar device, special attention was paid to the quality of a sample. In the study a variation 2 x 10^* in 
concentration over the entire volume was achieved which permitted scattering studies up to |t| > 10^"^. Within the 
interval 1.5 x 10"'^ < |r| < 10~^ the inverse correlation length was obtained from the width of the Lorentzian fits as 
a function of temperature, and the correlation length exponent was extracted from the power-law fits. The fit to the 
power ansats for x with background term yielded 7 after extrapolation to the wave vector length q = (see Table |) . 
An alternative substance with the base Mn (MnpZni_pF2) was studied by neutron scattering experiments a year later 
M. The sample with p = 0.75 had overall spread in concentration 0.001 and allowed to perform measurements up to 
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|r| ~ 4 X 10~^ (see figure ||), while the p = 0.5 sample was of lower quality with a spread in concentration of 0.005. 
Quite high quality of the samples as well as the temperature control up to 0.05 K allowed to obtain critical amplitudes 
ratios and exponents with good accuracy. However, systematic errors, which were assumed to be due to the resolution 
corrections, the quasi-elastic approximation, concentration fluctuations and background effects, appeared to be of 
more importance than statistical ones. The authors stressed that since no correction-to-scaling terms were used, the 
exponents found are effective ones. However the p — 0.75 sample was shown to exhibit RIM critical behaviour over 
the reduced temperature interval 4 x 10"'* < |r| < 2 x 10~^, whereas the critical exponents at higher dilution p — 0.5 
had not reached their asymptotic RIM values (see Table | and figure ^ . 

The results mentioned above for the single crystal of transition metal fiuorides were corroborated |^ by neutron 
scattering technique for a sample of dysprosium aluminum garnet Dy3Al5 0i2 powder. A powder was prepared of 
this cubic non-coUincar Ising antifcrromagnet together with non-magnetic Yttrium to avoid gradients in composition 
which unavoidably occur. The critical exponent /3 increased in disordered sample in comparison to the value of the 
"pure" Ising model (see Table ^ and provided an extra evidence of the mapping between experiment and theoretical 
insight. 

It is known that static critical behaviour of the diluted Ising magnet in a uniform field H along the uniaxial direction 
behaves according to the random-field Ising model p^ , |97t . The random-field Ising model is the subject of intensive 
recent experimental studies (see |9^,^ for a review). Two studies on the NS measurements in RIM appeared in the 
middle of 90th in the context of random-field model investigations ||l^,|l^ . In the NS studies the problem of the Bragg 
scattering saturation due to extinction effects was reduced by growing a 3.4 fim epitaxial film of Feo.5Zno.5F2 on a 
(001) ZnF2 substrate. The small X-ray rocking curve line-width of the (002) reflection showed the film to be of very 
high quality which was in part due to the nearly identical a-axis lattice parameters of the substrate and substance. 
Nevertheless, the film contained 10"* lattice spacings and could be considered as a three-dimensional object as was 
proven in FeF2 showing the critical behaviour of the "pure" d = 3 Ising magnet |99). The transverse (100) Bragg 
scattering scans data appeared to be well described by a Gaussian distribution with a background term. Fitting the 
Bragg amplitude vs temperature to the simple power law a critical exponent /3 — 0.35 was found for iJ = 0, however, 
a rounding of the expected critical behaviour of the magnetization close to Tn prevented an accurate analysis of the 
data. To reveal whether this was caused by poor sample quality the measurements on the 3.4 fim film of Feo.52Zno.48F2 
were performed ||l^ . The higher resolution data confirmed the value of f3 obtained and showed that rounding was only 
due to resolution. Recently in a thin crystal of 0.44 mm the scattering measurements were extended from |t| > 10~^ 
to |r| > 10"** to obtain the correlation length and susceptibility critical exponents [p^ , p^|j2C[| . 

The Mosbauer spectroscopy (MS) studies of the critical behaviour started in 1986 on a class of Fe-based substances 
with various magnetic atom concentrations H]. The method was tested earlier for the "pure" Ising antifcrromagnet 
FeF2, where the sublattice magnetization was shown to be proportional to the iron hyperfine field h which can be 
measured in MS [po| . To apply the approach to a dilute magnet, the variation of the concentration of Zn was reduced 
to 10~^ and a temperature stability of 0.002 K was achieved. In order to fix the critical temperature the values of 
P were chosen such that the plot h{Ty^^ vs T became a straight line intercepting the abscissa axis at T^r. This 
method appeared to be not sensitive to particular values of /? and gave the accuracy 0.05 K for critical temperature 
location. Though the critical exponent f3 can be obtained from the curve slope in double-logarithmic plot, this method 
produced results only after separating the data into two intervals in r which were characterized by obviously different 
critical exponents: one of the "pure" Ising model and the other one of RIM. The crossover in f3 occured within a 
very narrow range and at relatively large values of 10~^ > |t| > 10~^ for small dilution p > 0.93. The subsequent 
MS study appeared two years later ||T^ and is characterized by a very detailed analysis of the data. An advantage 
was the high quality of the sample grown from a stoichiometrical mixture of FeF2 and ZnF2 powders, each prepared 
by reacting a metal sponge with dry HF at -1-900°. In the experiment concentration gradients in a single crystal 
absorber were minimized by choosing the direction of the 7-ray parallel the growth axis, both perpendicular to the 
plane of the sample disk of exposing sizes 4 x 5 x 0.1 mm. 20 constant-acceleration spectra were obtained within 
3 X lO^'' < |r| < 0.86. Very good temperature stability 3 mK/day and account of correction-to-scaling permitted to 
obtain the asymptotic RIM magnetization critical exponent [3 from the data in the investigated interval (see Table |). 

The optical linear birefringence method (LB) introduced in Ref. |^ for the study of critical behaviour of the magnetic 
part of the specific heat Cm in disordered magnets was continued in Refs jl^,|lj,|l^ . The proportionality between the 
temperature derivative d{An)/dT of the optical birefringence An and Cm is assumed for studies of optically transparent 
materials. While the caloric properties in the critical region can be measured also by the pulsed heat technique, the 
LB method has several advantages: first of all, the non-magnetic contribution to the temperature derivative of the 
optical birefringence is insignificant, contrary to thermal techniques where the non-magnetic phonon background is 
often large and difficult to eliminate; secondly, one can minimize the concentration fiuctuations effects by applying a 
laser beam perpendicularly to the concentration gradients. In the first experiment jT^ a class of single crystals with 
various impurity concentration showed the exponent a to be independent of concentration and in very good agreement 
with the theoretical predictions (see Table III), while the two subsequent studies [|l^,|l5 showed the equivalence of 
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the LB method results w ith t he direct heat pulsing measurements. Magnetic X-ray scattering successfully applied to 
the "pure" Ising crystal |100| was applied to the "random" Ising crystal for high dilution in Ref. For samples of 
sizes 5.7 x 6.4 x 8.9 mm and of good crystallographic as well as mosaic quality the mean-field value for the exponent 
/3 has been found up to |t| < 0.06, while closer to Tjv an exponent (3 consistent with the RIM value was obtained. 



IV. MONTE CARLO SIMULATIONS 

The power-law singularities like in physical quantities only appear in the thermodynamic limit when the system 
volume and the number of particles tend to infinity. Therefore an obvious obstacles in computer "experiments" is 
that simulations can be done for a system of finite size only. Moreover the asymptotic temperature interval fails to be 
reached since Tc itself is not known exactly, whereas the location of the critical temperature is crucial for the accuracy 
of critical exponents determination. One can extrapolate to the thermodynamic limit assuming that the regime of a 
constant asymptotic critical exponent is established already starting from a certain finite system size and within an 
interval around Tc of a certain non-zero size. The last assumption exploited in computer simulations leads to a very 
narrow temperature intervals below and above Tc reliable for data sampling. The upper bounds are the temperatures 
where correction-to-scaling accounts for, while the lower bounds follow from finite-size effects and are of order 10"'^ 
for typical system sizes achieved now in the Monte Carlo (MC) simulations. Therefore, the exponents obtained in 
the manner described often are effective critical exponents which characterize the critical behaviour in the observable 
temperature range. On the other hand different time scales govern the dynamical behaviour of pure and disordered 
spin systems and only very long-running simulations yield reliable data. The relaxation time increases drastically not 
only when approaching Tc- Moreover, in the early simulations it increased also with dilution and yielded in several 
decades greater values when passing from a pure p — 1 system to a diluted one with p — 0.6. However an application 
of a more elaborated simulation technique resulted in an opposite behaviour: the decrease of the relaxation time for 



a certain disordered sample with an increase of p [102|. And it is the configurational averaging that leads to overall 
increase of the computation time. The statistical errors for thermodynamic observables result mainly from variance 
in configurational space and are larger than usual statistical errors in finite-size simulations with usual now statistics 
of 10^ MCS (Monte Carlo steps per spin). Consequently, the accuracy is bounded by the available CPU time from 
one side and limited by the number of samples to perform the configuration average of a disordered model [ [76[ from 
the other side. Efficient algorithms together with high-speed computers partially solve the task. 

Despite of this unpleasant situation Monte Carlo studies provided deep insight into the origin of the phase transition 
in the RIM as well as resulted in reliable numerical values for the critical exponents. MC studies of three dimensional 
RIM systems have been carried out over the last two decades . The first search |2l|] for universal critical char- 

acteristics of a RIM was performed on a simple cubic lattice of a size up to 30^ on the basis of an importance sampling 



MC method 1 101]. The positions of specific heat peaks corresponding to finite-lattice pseudo-critical temperature were 
extrapolated by means of finite-size scaling theory to Tc, and the accuracy of 500 — 5000 MCS per data point with 
averaging over several different starting configurations of impurities did not allow to discriminate between disorder 
relevance or irrelevance. From a fit the value of the order parameter exponent /3 and the susceptibility exponent 7 
of RIM could not be distinguished within the error bars from the corresponding exponents of the pure (non-diluted) 
system for any dilution (see Table II, Ref. |2l|). As a possible explanation it was stated that the lattice studied were 
too small to reach the critical region. The predictions of the authors that "impure" critical behaviour can be observed 
only for systems with larger critical regions (larger apure) was objected by the MC simulations of RIM on larger 
lattices with narrower dilution range [g^, higher statistics of 5000 11000 MCS, larger number of averaged samples 
and more detailed analysis of simulational data. The phase transition temperature Tc was located by adjusting data 
of specific-heat and numerical derivative of the energy. In this way the critical temperature accuracy was improved to 
0.004 and the relevance of disorder for the universal critical properties was concluded. This result was confirmed by 
magnetization data, though the critical exponent P was found to vary continuously with the magnetic sites concen- 
tration and no temperature region with constant P was found. Similarly to Ref. |2l| the authors of Ref. ||2^ stressed 
that while they cannot exclude the existence of a tiny impure critical region it must be unobservable for the weak 
dilution experimentally as well as by Monte Carlo methods, and effective critical exponents varying continuously with 
dilution would always show up according to the two fixed points scenario (see fig. Even the possibility of the 
existence of a line of fixed points, one for each impurity concentration, was not excluded (compare the dependence of 
the magnetization critical exponent f3 on the concentration of magnetic sites in Table |l| Ref. ||2^). 

The results mentioned above appeared to be in a good agreement with the MC data of Ref. |23|] obtained on the 
basis of a multi-spin coding program where twice as large systems and 8000 MCS statistics were used for a given 
concentration and temperature. The effective exponent /3 was found to increase continuously with dilution as well. 
However this was believed to be a consequence of the fact that one needs to compute the equilibrium magnetization 
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much closer to the critical point. Two years later the conclusion that reliable critical exponent values could not be 
obtained by the simulation of systems with sizes I < 20 was drawn in Ref . . Statistical errors in the determination 
of the critical temperature from susceptibility and/or specific heat data, together with the extr apola tion by finite-size 
scaling were thought to be responsible for this situation. In particular correction-to-scaling 103 should be taken 
into account. Instead, the simulation with the statistics 5000 11000 MCS, but no configuration average allowed 
a determination of the exponent /3 obviously larger than the pure one, while the determination of 7 were not done 
because of a large scatter of the data. A reanalysis of the previous MC simulations [pT|-p3t led the authors to the 
conclusion that a large part of critical regime is dominated by the crossover from the pure Ising exponents to RIM 
ones: the results of Ref. ||2l[] may correspond to a plateau value of the effective critical exponent governed by the pure 
Ising fixed point, while refined results of Refs [g2 23 1 correspond to an intermediate regime before the asymptotic 
regime, unaccessible near T^, is reached. Enforcing the concept of effective critical exponents ||6^,^ (sec formula 
(||)) the authors of Ref. |2^ claimed that the effective exponents they found were in qualitative agreement with the 
renormalization group flows (see fig. |^). 

Soon it became clear that the concentration dependence of the critical exponents observed previously in MC 
simulations may be also due to large relaxation time at criticality. The Swedsen-Wang algorithm applied to RIM 
simulation of Ref. resulted in the conclusion that a finite-size analysis of susceptibility data is inadequate to check 
the new critical exponents appearance for a RIM since the ratio j/v does not change with the system size. On the 
other hand, the determination of Tc from susceptibility maxima and the fourth-cumulant intersection method allowed 
to calculate the effective critical exponent 7 by fitting the data points in successive time intervals. In this way the 
susceptibility and the correlation len gth critical exponents were found to be independent of concentration in a wide 
range of dilution (see Table II, Ref. p^). As observed in Ref. |2^ the exponent 7 appeared to be higher than the 
theoretical asymptotic value. Thus it has been argued that Jeff is a non-monotonic function of reduced temperature 
with a maximum value, the larger the bigger the dilution. For the particular case of magnetic sites concentration 0.8 
the prediction was investigated in Ref. |2^ by the Swedsen-Wang algorithm and the single-cluster generalization of 
Wolf for a rather large system. While the j/iy ratio confirmed the result of Ref. |2^], the effective critical exponent 7 
did not show a maximum but increased continuously approaching Tc- 

The assumption that the interaction among block-spins in diluted systems may be represented by renormalized 
couplings of pure systems was used in the MC renormalization group approach for RIM in the weak dilution region 
(p > 0.8) The values for the correlation length critical exponent v were found to depend on the concentrations 
of the magnetic sites. For p = 0.8 ly satisfied the exact inequality v > 2/d for a diluted system |Q and this led 
the authors of Ref. to the assumption that they obtained the asymptotic value of the correlation length critical 
exponent. They subsequently estimated the width of the asymptotic critical region for this dilution (Irjcrit < 3 • 10"'*) 
where no influence of the pure Ising model fixed point is seen. Alternatively, for p — 0.9 the inequality [72) is not 
satisfied for v and either the critical region Irjcrit < 1-3 ■ 10"'* was not reached and thus the critical behaviour is that 
of pure Ising model, or the crossover to the RIM asymptotical critical behaviour for p = 0.9 is observable only in 
larger systems. For p = 0.4 the accuracy of data obtained was too poor to find any values. 

Precise MC data of Ref. ||2^ became availa ble by application of the new cluster algorithm 104 and a refined 
vectorized implementation of local algorithnisj 105 1 . A progress in MC studies occurred when an improved version 
of the multi-spin coding program appeared | 105^ |. It allowed to verify all previously obtained simulation data for 
the RIM with the unprecedented statistics of up to 3 x 10^ MCS and 10 averaged configurations for systems with 
p > 0.8 and up to 1.2 x 10^ MCS for higher dilution. For instance, in the test simulation the statistics allowed to 
obtain the susceptibility and magnetization critical exponents of a 2d pure Ising model with an accuracy of 1% and 
3% respectively. For the RIM, the effective exponents 7,/? and ( — I — P (the last one describes the divergence of 



29| and were shown to be concentration dependent in 
ue as well as 7 were obtained by fitting susceptibility 



the magnetization-energy correlation function) were obtained 
the concentration region 0.5 < p < 1 (see Table |l|). The Tc va 
data to a simple power law; the Tc location was then verified by the 4th and 6th cumulant of magnetization. The 
investigation of the energy-magnetization correlation function allowed to check if scaling holds; the critical exponents 
C and P were obtained by fitting the data to the simple power law with the value Tc taken from the susceptibility 
analysis. All data showed power-law behaviour within the chosen temperature range but the constant values of critical 
exponents changed with concentration (see figure ^). For instance, 7 increased from the pure Ising value and smoothly 
changed with dilution achieving plateau value at p — 0.5, ^similarly increases while ^ decreases in such a way that 
their sum equals 1 within errors. The conclusion of Ref. ||29| states that while new critical exponents changing with 
the dilution were observed, no line of stable fixed points exists as supposed e.g. in The result originates from a 
crossover from pure to diluted and percolation to diluted regimes with effective exponents for all concentrations. 

The evidence of crossover phenomena motivated the authors to undertake a more systematic study. The data of 
Ref. were revised three years later in Ref. |30| where a cumulant method was used in order to determine Tc- For 
thermal averaging up to 40000 MCS were performed and for the configurational averaging up to 32 configurations 
were used (compare this with 5000 MCS and only several samples for the configurational average in Ref. |2l)). The 
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critical exponents C, P and 7 were shown to vary with concentration in consistency with the two fixed point scenario 
for the weak dilution. However for the strong dilution (p = 0.5, 0.6) an influence of another "percolative" fixed point 
was assumed. A central part of Ref. |^ dealt with the magnetization-energy scaling function. It was shown that 
dilute systems have a complex crossover behaviour before they reach their asymptotic critical region with values 
for the exponents consistent with the weak random fixed point. For instance, pure systems reach their asymptotic 
limit for small system size. In weakly random systcms(p > 0.8) the asymptotic values for the exponents are smoothly 
approached from below whereas in strongly diluted systems from above. This happens at the characteristic system size 
Xp depending on the concentration p. Xp was estimated as 20 30 lattice constants for pure, 50 for p — 0.95, 0.9 and 
about 100-^ 150 for p — 0.6; at the percolation point Xp diverges. Thus more strongly disordered systems necessitate 
a refined analysis with an appropriate treatment of their percolative structural effects relevant for non-asymptotic 
sizes and correlation lengths. Accidentally the crossover function changes its sign at p ^ 0.8 so the asymptotic values 
appear to be reached already for small systems (c. f. Refs p7| , p9t ). 

Due to works |^,^ and especially [ p9pC| ] it became clear that the concentration dependent critical exponents 
found in MC simulations are effective ones, characterizing the approach to the asymptotic region. This point of view 
found its support in a conjecture about a step- like universality of the 3d diluted magnets An attempt to study 
the diluted model at p — 0.6 by means of taking MC data from different subsystems of one large system was made 
in Ref. |^ . It was concluded that the resulting values of the critical exponents are strongly influenced by the size of 
the system. 

Recently, the critical behaviour of the RIM was reexamined by MC method in Refs [M|J3^] for concentrations 
p = 0.8,0.6 and in Ref. |Q for 0.4 < p < 0.9 respectively. In particular, the simulations |35 revealed that a disorder 
realized in a canonical manner (fixing the fraction of magnetic sites) leads to different results in comparison with 
disorder realized in a grand-canonical manner (see Table 0). Studies of Ref. [|| were based on the importance of 
taking into account the leading correction-to-scaling term in the infinite volume extrapolation of the MC data and 
thus the analysis does not agree with the data of Ref. The results of the simulations for the concentrations 

p = 0.9, 0.8, 0.6, 0.4 were extrapolated to infinite system size (see figure 0) and lead to the proof of universal critical 
behaviour of the site-diluted Ising model in a wide range of concentrations. In particular the value of the correction- 
to-scaling exponent to was found to be a; = 0.37 ± 0.06 which is almost half as large as the corresponding value in 
the pure d = 3 Ising model ||7^ uj = 0.799 ± 0.011. The smaller the value of uj the larger the interval where it has 
to be taken into account (c.f. formula (|^) of the chapter ||). Thus the smallness of uj in the dilute case explains its 
importance for an analysis of the asymptotic critical behaviour. 

Another important question considered in Refs |^,^ is the problem of self-averaging in RIM. The Gibbs approach 
to the static collective phenomena rests on the statistical independence of macrosamples according to the short range 
nature of inter-particle interactions. In agreement with this approach any thermodynamic extensive quantity M is 
(strongly) self-averaged. It means that the normalized square width Rm of the squared variance of its sub-sample 
values behaves as Rm = Vm/M^ ^ 1/n ^ where n is the number of sub-samples and I is the system linear size. 
However, in the vicinity of the critical point the statistical independence does not hold since the correlation length of 
a system ^ can be arbitrary large ^ ~ and thus subsystems cannot be considered as independent. The concept of 
weak self-averaging takes account of the case when a number xi,0 < xi < d exists such that Rm at criticality scales 
as In contrary, if Rm const 7^ 0, M is called to be not self-averaging. It has been predicted on the basis of 

heuristic arguments that for random models all extensive quantities are strong self-averaging far from criticality. Yet 
for a quantity M scaling at the critical point as F the strong self-averaging should fail. Here the squared variance 
Vm is expected to scale as /^p+a/f f-Qj. arandom < implying Rm l"^" or weak self-averaging |10(:]. 

Lack of self-averaging in RIM is not only of high theoretical interest. The reliability of M C simulations depends 
on the answer whether the increase of the lattice size improves the statistics of the simulations. If a quantity is not 
self-averaging the simulational data are unreliable. Theoretical studies based on the renormalization group approach 
confirmed the strong self-averaging for / ^ ^. In contrast, special MC investigations in finite systems found no 
self-averaging in th e ca se of relevant disorder a^^^^^ > , while weak self-averaging appeared to be the case only for 
irrelevant disorder [107|, in disagreement with Ref. |l06|] . The MC simulations of the Refs p^ , ^ were performed in 
order to resolve this puzzle. It was shown that the normalized square width Rm goes to a constant for large I, which 
is independent from the dilution of grand-canonical type, while this is not the case for the canonical type |34[| . The 
last result, however, may have its explanation in the very slow approach of R to its universal asymptotic value in the 
case of canonical realization of 3d random Ising model, estimated as l"^'^ p6| . 

The evolution of the self-averaging from pure Ising model to RIM has been studied recently in order to determine 
the transition zone between pure and diluted Ising models' universality classes. It was shown that the transition zone is 
smoothly dependent on concentration of magnetic sites and independent of lattice size. In contrary, critical exponents 
values did not depend on concentration and were found identical to previous data of Ref. [3^ ] . The universal value of 
the normalized square width of the susceptibility in the infinite-volume limit was estimated to equal i?^(oo) = 0.155. 

Apart from short-range site-dilution other realizations of disorder have become a subject for MC simulations recently. 
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Thermally diluted Ising model has been studied as a generalization of the RIM There, the realization of the 

vacancy distribution is determined from the local distribution of spins in a pure Ising model at criticality. The critical 
properties, in particular, the universality class of the model appeared to differ strongly from the RIM one, however, 
they agreed well with the theoretical predictions for long-range-correlated disorder [ 

The problem whether the RIM fixed point also describes the phase transition in the Ising model with random 
bonds has been studied explicitly in Refs [^0. Using a numerical renormalization-group analysis the RG flows for 
random Ising models have been obtained and the existence of a fixed point characterizing the random Ising model 
irrespectively to the type of disorder has been shown . 



V. RENORMALIZATION GROUP THEORY EXPANSIONS 



In the subsequent sections and VI we will review results on RIM critical exponents obtained my means of 
renormalization group (RG) methods. In section ^ we report the main relations of the field theoretical RG approach 
and dwell upon perturbation expansion series available. Section VI deals with the RG series resummation methods 
and the results obtained on their basis. 



A. Renormalization 



To describe theoretically the long-distance properties arising in different systems in the vi cinity of the second order 
phase transition point it is standard now to use a field theoretical RG approach [111 112]. The renormalization is 



used to remove divergences that occur during evaluation of the bare vertex functions in the asymptotic limit. For the 
RIM one-particle irreducible bare vertex functions are defined as: 

S{qi + ■■■ + qN)r^{qu qN^m^, Uo, Vol Aq) = J e'^'^'^' + -+'^-^-\^{n) ■ ■ ■ ^(rw))^p-™dri • • • dr^, (12) 

where the angular brackets denote statistical average over the Gibbs distribution with the Hamiltonian (H) in the 
replica limit n — > and the subscript |ipi indicates that only one-particle irreducible diagrams are taken into account. 
The functions depend on sets of momenta qi,...,qN (with Ag as a momentum cutoff) and the bare parameters 
r7io,uoi^^o of the Hamiltonian (H). Divergences in (|l^ ) occur in the asymptotic limit Aq — *■ c». Their removing 
is achieved by a controlled rearrangement of the series for the vertex functions. Several asymptotically equivalent 
procedures serve this purpose. For the purpose of this presentation we will use two complementary approaches: (i) 
dimensional regularization and the minimal subtraction schem e Ill3[ and the fixed dimension renormalization at zero 
external momenta and non-zero mass (a massive RG scheme) |114| . 

Let us formulate relations of the renormalized theory. The renormalized fields, mass and couplings (j),m,u,v are 
introduced by: 

^ = (13) 

"^0 ^ ^m^m^ (14) 



Here, e = 4 — rf, /i is a scale parameter equa l to the renormalized mass at which the massive scheme |114| is evaluated 



or in the minimal subtraction scheme |113] it sets the scale of the external momenta. Z^, ^m^, Z/^^u, ^4,t; are the 
renormalizing factors. The renormalized vertex functions F;^ expressed in terms of the bare vertex function by: 

F^(gi, . . . ,gAr;m,u,w) = Z^^^Tq (gi, . . . , ^at; mo, Mq, "fo)- (17) 

are finite. This is the main content of the multiplicative renormalizability of the field theory defined by the Hamiltonian 

(i-. . . ^ ^ . 

First let us consider the minimal subtraction scheme. Here, the renormalizing Z-factors (Uy)-dlq) are determined 
by the condition that all poles at e = are removed from the renormalized vertex functions. The RG equations are 
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written bearing in mind that the bare vertex functions do not depend on the scale and therefore their 

derivatives with respect to fj, at fixed bare parameters are equal to zero. So one gets 

/^|;r-|o^^|-z;-/^r-|o = o, (is) 

where the index |o means a differentiation at fixed bare parameters. Then the RG equation for the rcnormalized 
vertex function reads: 

d d d d N \ 

tJ-TT + Pu-^ + Pv-^ +"/rnm- -"f^)T^{m,u,v,fi) = 0, (19) 

ou ov dm 2 J 

and the RG functions are given by 

(20) 
(21) 
(22) 

1 d\nZ-\, , , 

Using a method of characteristics 115 to solve the partial differential equation ( p^ ) one can write a formal solution 
as: 

r^(TO, V, fi) = X{e f/^r^{Y{£)m, u{i),v{£), fi£), (24) 
where the characteristics are the solutions of the ordinary differential equations (flow equations) : 

tj^u{l) = P^{u{l),v{£)), l±v{l)=p,{u{t),v{l)) (25) 

with 

X(l) = r(l) = l, u(l) = M, v{l)^v. (26) 

For small values off, equation ( p^ ) maps the large length scales (the critical region) to the non-critical point £ — I. 
In this limit the scale-dependent values of the couplings u{£), v{£) will approach the stable fixed point, provided such 
a fixed point exists. The fixed points m*, v* of the differential equations ( p5[ ) are given by the solutions of the system 
of equations: 

/3„(u*, V*) = 0, 

l3y{u\v*)^Q. (27) 
The stable fixed point is defined as the fixed point where the stability matrix 

Bij , Ui = {u, v} (28) 

possesses eigenvalues uji,uj2 with positive real parts. The stable fixed point, which is reached starting from the initial 
values ^0 in the limit £ ^ corresponds to the critical point of the system. In the limit f — > (corresponding to 
the limit of an infinite correlation length) the rcnormalized couplings reach their fixed point values and the critical 
exponents rj and v are then given by 

77 = 70(^*,«*), (29) 
llv = 2{l~^.m{u\v*)). (30) 
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In the non-asymptotic region but near the fixed point deviations from the power laws with the fixed point values of 
the critical exponents are governed by the correction-to-scaling exponent 



uj = min{uJi,LU2) (31) 
in accor danc e with the Wegner expansion (||). The rest of the critical exponents are obtained by familiar scaling 



laws |67,111|: 



a = 2-di/, /3= ^(d-2 + 77), 7 = 1/(2-77) (32) 

which can be shown to hold from the solutions (|2^). 

The flow equations ( psf ) can serve to describe the approach to criticality in a larger region where corrections to scaling 
do not suffice. As it was mentioned in chapter | out of the asymptotic region physical observables are characterized by 
effective exponents introduced to describe a crossover from the background behaviour to the asymptotic critical one. 
In the RG language they depend on the flow parameter £{t) through the dependence of couplings on £. In particular, 
according to the definition (see formula (^) for the magnetic susceptibility effective exponent 7e// we have: 

(r) = = liu{iiT)),v{iir))) + ..., (33) 

where the second part is proportional to the /3-functions and comes from the change of the amplitude part of the 
susceptibility. This part is natural to neglect in the vicin ity of a fixed point. Moreover, the contribution of the 



amplitude function to the crossover does seem to be small pi6| , |ll7| . Under this restriction, the effective exponents 
are simply given by the expressions for the asymptotic exponents (pG) but with replacing the fixed point values of the 
couplings u*,v* by the solutions of the flow equations (|25|): 

l/,,^„(i)^2{l~jMe),vm- (34) 

In the massive RG scheme the Z-factors (p^-(p^ are calculated from the vertex functions ( p^ ) at zero external 
momenta gi, . . . , gjv and non-zero mass at fixed space dimension in the limit Ao — > 00 . This normalization conditions 
lead to the equation for the rcnormalized vertex functions which is called Callan-Symanzik equation. Differentiation 
of the proper bare vertex function Tq with respect to the rcnormalized mass (c.f. eq. (p^) gives: 

again, the index |o means a differentiation at fixed bare parameters and a new vertex function 

appears. This vertex function differs from Tq by an extra term ^p'^iR) inside the averaging < . . . > in (p^. As a result, 
for the renormalized vertex functions one gets an inhomogeneous Callan-Symanzik equation containing F^'^' in 
the right-hand side. However, close to the critical point to = the right-hand side can be neglected with respect to 
the left-hand side and one arrives to the homogeneous Callan-Symanzik equation which repeats the structure of the 
RG equation (^: 

where the coefficients 7^ are defined by the relations ([20|)-(|2^), but the parameter ^ there is to be understood 

as the renormalized mass to. The partial differential equation { p7[ ) is solved again by the characteristics method 
sketched by the relations (^4|)-(j2^) and leads to the fixed point relations given by (|2^)-(|3l|). 

For the sake of completeness let us note that the finiteness of the renormalized vertex function with one 0^ insertion 

(12)'' . ■ . . ~ 

Tj^ ' is achieved by the familiar renormalizing factor Z^2 : 

Tfl^'^^ {k; q, -q; to, m, v) = Z^2r'^^'^^ (fc, q, -q, toq, wq, ^o)- (38) 
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Then, formula 

27m = 70+70=^: 



^) for the correlation length critical exponent v may be recast in terms of Z^i by a substitution 



which follows from the relations = Za,iZ^, and 



702 



fid \n Z ^2 / d fi\o ■ This leads to the relation: 

(39) 



The explicit expressions for the RG (3- and 7- functions are scheme dependent: they do differ in different renormal- 
ization schemes. Subsequently the fixed point coordinates are scheme dependent as well. However, the RG functions 
coincide in different schemes once calculated at the fixed point. This leads to the same a symptotic values of the 
critical exponents: they are universal and do not depend on renormalization scheme 111 112]. In the next subsection 
we will provide available expansions for the RG functions of the RIM. 



B. Perturbation expansion series and their "naive" analysis 



Expressions for the RG functions of the RIM are obtained as series in the renormalized couplings u and v. The 
perturbation theory in powers of m, v is in fact a perturbation theory in number of integrations in fc-space which, 
on the other side, cor resp onds to the number of loops in Feynman diagrams in the diagrammatic representation for 
the vertex functions Ullll . By now, the RIM RG functions are known up to the order of 5 loops in the minimal 
subtraction scheme |118|] and with record 6 loop accuracy when calculated directly for d = 3 in the massive scheme 

. . . . 

Written in the minimal subtraction scheme the functions read: 

/3„ = -u (e _ u - 3/2 V + 17/27 + 23/12 uv + 41/32 + . . . + f3^^^^^), (40) 

Py = -^;(£-^;-2/3^i + 21/32^;2 + ll/12^>M + 5/27^i^ (41) 

70 ^ 1/54 + 1/24 + 1/64 + ... + 7^^^^^ , (42) 

702 = l/3u+ l/4w- l/9u2_ 1/4 3/32^2 _^ _ ^ ^i^LA) _ ^^^^ 

Here, we write down the functions only up to the two-loop approximation (2LA), referring to the paper | |5^ where 
they were obtained in 3LA. Four- and five-loop contributions may be derived from the RG functions of the anisotropic 
cubic model obtained in Ref. |118{ . 

Note, that in the minimal subtraction scheme the dependence on the space dimension d is trivial and enters 
expressions (|40|)-(^3|) only via one single term proportional to e = 4— d explicitly written in the /?- functions (|40|), (|4l]). 
On contrary, in the massive scheme the space dimension d enters the expressions for the loop integrals, corresponding 
to each Feyn man diagram of perturbation expansion. Subsequently the theory is evaluated at space dimension of 



interest [121|. The series for RG functions of the d — 3 RIM read: 

/3„ = -u (1 - M - 3/2 w + 308/729 + 104/81 uv + 185/216 + . . . + , (44) 

py=-v{l-v-2/3u + 95/216 + 50/81 vu + 92/792 + . . . + Pi^^^'^) , (45) 

70 = 8/729 + 2/81 uv + 1/108 v^ + . . . + 7^^^^^ , (46) 

702 = l/3w+ l/4w- 2/27^2- l/6ut;- 1/16^2 + . . . + 7^^^"^^ . (47) 

Again, we display here only the two-loop functions. In the e expansion these functions were first obtained in Ref. 
psf . Three-loop terms first reported in Ref. [|8| contained some errors, partially corrected in Ref. Finaly free 
of errors three-loop expressions were reported in Ref. fQ. Subsequently four-loop series were obtained in Ref. p^ ] 
and only recently five-loop [p3 and record six-loop [|66|| expansi ons became available. Universal critical amplitude 



ratios at d = 3 were first obtained in three-loop approximation | 119 | and are known by now within the accuracy of 



five loops 1120 1 



As we noted above, the massive RG scheme does not necessarily mean evaluation at d = 3. A possibility to apply 
the scheme in order to get the RIM RG functions at arbitrary non-integer space dimension was outlined in Ref. . 
For the two-loop RG functions one gets Q: 

Pu{u,v) - -(4 - d)u[l -u~^v + -^ [9(ii - i) + + ^ [I2{h - ^) + i2]uv + i [21{n - i) + ia]^'}, (48) 
/3„(u, v) = -{4- d)«|l - w - \ + i [ll(ii - 1) + i^] v^ + ^ [6(zi - i) + 12] «w + ^ [3(ii ~ i) + za] M^}, 
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702 (u,w) 



= -2(4-rf){ 
= (4-d){i^ 



27 
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1 



16 
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(49) 



27 



12 



32 



The space dimension d enters expressions (^8|)-(|50|) also by the rf-d epen dent two-loop integrals ii{d), 12 (d). Their 
dependence on d is shown in the figure ||. Evaluating the integrals |121] for d = 3: ji(3) = 2/3, 12(3) = —2/27 one 
recovers the two-loop contributions of the expressions (|4^)-(|47|). In next order of perturbation theory f or no n-integer 
d the RIM was derived in Rcf. |6^, the values of the corresponding loop integrals are obtained in Ref. [122|. 

With the series for the RG functions at hand there are two different ways to proceed in any of the two RG 
approaches (p0|)-(p3|) or (Q)-(^). Subsequently, this leads to four differe nt sc hemes for our analysis. Indeed , th e 
minimal subtraction scheme can be realized by the familiar e-expansion |123| as well as directly at d = 3 [124|. 



Similarly, the expressions for the RG funct ions obtained in the massive scheme can be directly solved for the fixed 
points (FP) or by the pseudo-e expansion |l25|| . However, as we will see below, only two of these schemes lead to 
reliable results for the RIM critical behaviour. Before we proceed let us discuss this situation in more details. 

• 1. First let us consider the functions (|40| )-( ] 43| ) and perform the e-expansion. In order to obtain the expansion, 
one should (i) solve the fixed point equations (|27|)with e as a small parameter and then (ii) substitute the coordinates 
u*{e), v*{e) as series in e into the RG-7-function for the critical exponents for critical exponents. 

Looking for the fixed point solutions in one-loop approximation (i.e. leaving only linear in u, v terms in the brackets 
in (pO|), (^)) one finds three fixed point (compare with the fig. ||): the Gaussian fixed point G u* = v* = 0, the pure 
Ising fixed point I w* ^ 0, u* = (it describes n non-interacting Ising models: c.f. Haniiltonian (H) with v = 0) and 
the "polymer" fixed point P m* — 0, w* 7^ (putting m* = in (^) one gets the 0{n = 0)-model describing the scaling 
properties of self-avoiding walks). It is straightforward to check that the fixed points G and I are unstable whereas 
the fixed point P is stable but as far as u* > it is unaccessible for the initial values of couplings of our model. One 
does not encounter in ILA the fixed point R with both non-zero coordinates u* 7^ 0, w* 7^ 0: this is because the 
system of equations for the fixed points is degenerate on the one-loop level . This fixed point appears in the 

next, two-loop approximation and leads to the qualitative picture shown in the figure ^j. However, the degeneracy of 
the one-loop equations leads to the result, that instead of expanding in e one has to expand in y/e [^3|4^]. Proceeding 
as usual this leads to a y/e expansions for the critical exponents and the stability matrix (pM) eigenvalues |E^,El| : 



iy = 0.5 + 0.084115826^/2 - 0.01663203e -I- 0.047753516^/^ + 0.27258431e2 , 
?7 = -0.00943396£ + 0.03494350£3/2 _ o.04486498£2 -I- 0.02157321e^/2 , 
7 1 + 0.16823164£i/2 - 0.02854708£ + 0.0788288l£^/2 + 0.56450490£2 , 
coi = 2£ + 3. 704011194£3/2 -I- 11.30873837 £^ 

UJ2 = 0.6729265850 e^/^ ~ 1.925509085 e - 0.5725251806 e^^^ ~ 13.93125952 . 



(50) 
(51) 
(52) 
(53) 
(54) 



Two-loop expressions for the exponents were obtained in Ref. [ p5[ , three-loop results were presented independently 
in Refs |46[ and ^7\ . The -y£-expansion series for the RIM amplitudes ratio s [|127| is also available in the three-loop 
accurac3r [[l28[| . Due to the five-loop results for an anisotropic cubic model |ll§f| it was possible to obtain four- and 
five-loop v^-expansion for the critical exponents (p0|)-(p2|) and the stability matrix eigenvalues |6|] (|5^)-(|5^) of 
the RIM. 

• 2. T he se cond way of calculation which may be used in the minimal subtraction RG approach is the d — 3 
technique |124]. It consists in (i) fixing the value of £ = 1 in (40), (^), (ii) solving the system of fixed point equations 
( p7| ) numerically and (iii) substituting the numerical values of the fixed point coordinates into the series for the critical 
exponents. The RIM /3-functions are shown in figure ^ in two-loop approximation. One can see that they do not 
possess non-trivial fixed points at all! However, such behaviour is not surprising and is not a particular feature of the 
RIM. It is well known that perturbation series for the RG functions in the weak coupling limit are asymptotic at best 
and should be evaluated by means of special resummation technique. This will be the subject of the next Section 
VI, and here, to complete the list of possible calculation schemes we mention two more schemes in the massive RG 
approach: 

• 3. The massive RG scheme is implemented by (i) solving the system of fixed point equations ( |27| ) numerically 
and then (ii) substituting the fixed point coordinates numerical values into the series for the critical exponents. This 
is the way which most of the papers are devoted to in the RIM RG an alysi s (see the next section for more details). 

• 4. The massive RG scheme implemented by a pseudo-£ expansion [125|. This scheme is based on the observation 
that in order to analyze the series (44)-(|4^) one can imitate an £-expansion introducing an auxiliary parameter r 
instead of the zeroth-order term 1 in the massive scheme /3-functions (0), (^^. Then FPs and critical exponents 
are obtained as series in t with substitution t = 1 to obtain their final values. The described pseudo-£ expansion 
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|125| allows to avoid cumulation of the errors for the critical exponents from the error s of t he fixed point coordinates 
and the 7- functions. It proved to be highly efficient for the pure d = 3 Ising model |12£,7^. For the same reasons 
as the e-expansion for RIM turns to the ^/e expansion, the expansion in r turns to a -yr-expansion. On the basis of 
the six-loop expansions (^)-(^7|) we get the following y^-expansion for the RIM critical exponents and the stability 
matrix eigenvalues: 



1/ = 0.5 + 0.10291260 r^/^ + 0.01251853 r + 0.01270178 r^/^ + 0.05663757 + 0.03694322 r^/^ ^ 

77 = -0.00836820 t + 0.02173733 t^/^ _ 0.01487714 + 0.01733771 r^/^ , 

7 = 1 + 0.20582521 t^/^ o.02922117t + 0.01539608 t^/^ + 0.11858141 + 0.06658280 r^/^ , 



ui = 2t + 2.59761132t 



3/2 



7.51800557 + 39.86825804 r 



5/2 



(55) 
(56) 
(57) 
(58) 



UJ2 = 0.82330084 1.74713206 r- 1.26569350 r^/^ - 8.75074159 

These expressions are to be compared with the formulas (^-(|4|). 



40.98838378 T 



5/2 



VI. SERIES RESUMMATION AND NUMERICAL RESULTS 



A. Summability of RIM divergent series 



The core of the perturbation approach described above is the account of higher perturbation contributions. However, 
such a way does not lead necessarily to cumulation of the calculation accuracy by a "naive" summation of successive 
perturbation theory terms. It is well known by now that the weak coupling expansion series for RG functions have 
zero radii of convergence and are asymptotic at best [112|. Appropriate resummation procedures are to be applied 
in order to handle them. Early studies of critical phenomena by the £-expansion technique lead to a concept of a 
series "optimal truncation" as a maximal number of terms possessing convergent behaviour. Such behaviour is typical 
for asymptotic series [13C| where series expansion coefficients grow factorially. For example, the expansion of a RG 
function f{u) of pure Ising model in powers of a (single) renormalized coupling u: 



kU 



(59) 



was shown to possess the following behaviour 



A, 



ck'">{-a)''k\[l + 0(l/k)], 



(60) 



with known values 1 132 -134] of constants c, bo, a. The property ( |60| ) indicates the Borel summability of the series 
( p9| ) . The Borel resummation procedure takes into account the asymptotic behaviour of the coefficients ( |60| ) and maps 
the asymptotic series ( ^9|) to a conv ergen t one with the same asymptotic limit. We will describe the procedure in 



more details in subsections VI B and VI C . 

It is worthwhile to note, that there does not exists a proof of Borel-summability for the e-expansion of a pure 
Ising model. Only expansion in coupling (|5^) is proved to posses this property [131|. Nonetheless the rich amount of 
numerical estimates obtained on the basis of resummed e-expansion for the pure Ising model (see e.g. [112|) convinces 
one in its reliability. However it is not the case for the -ye-expansion of RIM. The fact that e-expansion will not be a ble 
to give information on critical exponents in system with quenched disorder was predicted already in Refs | |135| ,136| by 
studying the randomly diluted model in zero dimensions. The perturbation theory series of this toy-model appeared 
to be Borel non-summable. Moreover, such properties were shown |136] to be the direct consequence of the existence 
of Griffiths-like singularities [|o| caused by zeroes of the partition function of thepure system. Although the ^/e- 
expansion allowed to predict qualitatively new critical behaviour of the RIM |^,^ it seems to be of no use for a 
quantitative analysis. Moreover, from naively adding the successive perturbational contributions in the y^-expansion 
for the RIM stability matrix eigenvalues (|5^), (54) one observes that already in the three-loop approximation (~ e) 
UJ2 becomes negative and therefore no stable fixed point exists in strict y/e-expansion |^! Even the resummation 
procedures applied do not change this picture . Quite different convergence properties of the expansions for pure 
d = 3 Ising model and RIM may be seen already applying simple Fade analysis as it was shown in Refs ||59| , |63|| . 

The above mentioned divergent properties of the -y/e-expansion concern also the pseudo-y^ expansion (|55|)-([5^) 
derived in Section 0. So there remain only two out of four different ways of numerical analysis of the RIM RG 
functions. They are denoted by • 2 and • 3 in the preceding subsection |VB| . Both are based on analysis of RG series 
in two couplings. The nature of the RIM RG functions series in couplings u, v remains opened. Nonetheless, the 
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resummation procedures which in different modifications are used in analysis of asymptotic series have been apphed 
fruitfully to the RIM RG series as well |]5l|- |57| , |59| -|64| . The method which prevails in the resummation of the RIM RG 
series is the generalization of a Pade-Borel resummation t echniq ue |137| for the case of two variables. We will review 



results based on such resummation in the next subsection VI B 



Recently, it was shown that the expansions for the quenc hed diluted Ising model in = dimensions are Borel 
summable, provided a specific way of summation is applied [138|. This allowed authors of Rcf. both to recover 
the asymptotic behaviour of expansion coefficient as well as to apply to RIM six-loop massive RG functions ||6q| th e 
resummation technique refined by a conformal mapping. This method will be reviewed in the last subsection ]VI C|. 



B. Pade-Borel oriented resummation 

In order to perform the Pade-Borel resummation of a truncated (asymptotic) series in a single variable (e. g. series 
( p9| ) with L terms available) one defines the Borel-Leroy image of the initial sum S = J2i=o ^j'"* by: 

z— 

where r(a;) is the Euler's gamma function and q is an arbitrary non-negative number, which will serve as a fit 
parameter. The result ( |6l| ) is extrapolated by means of a rational approximant [M/N] (m), which is the quotient of 
two polynomials in u with M as the order of the numerator and N as that of the denominator (a Pade approximant 



[139|). Subsequently, the resummed function 5™'' is obtained in the form: 

POO 

S'^'iu)^ dtexpi-ty^M/Njiut). (62) 



In the single- variable RG functions analysis this procedure was initiated in Refs 137 



There are different ways to generalize single variable Pade app roximants for the case of two variables (u, u) |13!;] 



and, subsequently to generalize the resummation procedure (61), (|62|). One possibility is to construct the series in 



a single dummy variable x (so-called resolvent series [140 ). The variable is introduced by a substitution u —^ ux, 



V —>■ vx and x must be put equal to 1 to obtain final results. In accordance, the Borel-Leroy image of a truncated 
series S — J2o<i+j<L '^i.j'^'"^'' defined by: 



Then the series in x is resummed by formula (|62|) and evaluated at x = 1. We will refer to this method as to the 
Pade-Borel method similarly to the one-coupling case. 



Another way to proceed is to make use of the Chisholm approximants 1 141. 139] which are the generalization of 
Pade-approximants to many-variable case. A Chisholm approximant can be defined as a ratio of two polynomials 
both in variables u and v, of degree M and N such that the first terms of its expansion are equal to those of the 
function which is approximated. Again, the resummation is performed by (|6^ ) with Chisholm approximant instead 
the Pade one. This method will be referred below as the Chisholm-Borel resummation. 

For a given Borel-Leroy image many approximants, both Pade and Chisholm, can be constructed. However, 
restrictions are imposed naturally. Firstly, an approximant should be chosen in the form re-constituting the sign- 
alternating high-order behaviour of the general term of /3- and 7-functions [^. The approximant generating a 
sign-alternating series might be chosen in a form [Af/1] with the positive coefficients at the va riabl e x (or u and v) 



However, these are the diagonal approximants which give the best convergence in Pade analysis |139| (see also |63 for 
a toy model example). On the other hand, high degree of a polynomial in the denominator often leads to poles on the 
positive semi-axis. One can use analytic continuation and calculate the princ ipal value of the integral (p2|) to process 



the singularities however practical calculations reject the generalization |137|. The reason is both unexpected shift of 
fixed point location and different topological structure of the lines of zeros for the resummed /3-functions. The choice of 
a Chisholm approximant is even more complicated since often its coefficients are undetermined. For instance, already 
two-variable polynomials of a second order require two additional conditions to be imposed. Normally they are chosen 
to preserve certain symmetry properties. Provided the couplings uq and vq enter the Hamiltonian (H) symmetrically, 
the approximants must be symmetric in variables u and v in order not to introduce an additional symmetry. Another 
point is that by substitution t; = all the equations which describe the critical behaviour of the diluted model 
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are converted into appropriate equations of the pure model. However, if pure model is solved independently, the 
resummation technique with the application of Pade approximant is used. Thus, Chisholm approximant is to be 
chosen in such a way that, by putting any of u or w equal to zero, one reproduces familiar results for the one- variable 
case. This also implies a special choice of additional conditions. 

Most numerical results on the universal characteristics of RIM at criticality were obtained on the basis of the massive 
renormalization scheme by numerically solving fixed point equations (|27| ) for the rcsummed /3- functions (Q), (|4^) and 
resumming the 7-functions (|4^),(^7|) in the stable fixed point (scheme • 3 in the subsection VB). The study of the 
massive /3-functions of RIM resummed in this way revealed that starting from the two-loop approximation the random 
fixed point (R: see figure ^ is stable and is present in all orders of perturbation theory |51-56,6^, 32,34 142 1. As an 
example, in figures 10, ^we show the lines of zeros of functions I3y{u,v) in different orders of perturbation 

theory without resummation and resummed by the Chisholm-Borel method. One can see that without resummation 
all non-trivial fixed points are obtained only in the three-loop approximation. Resummation restores presence of 
non-trivial fixed points. 

The massive scheme of field-theoretical renormalization group was the basis of first numerical estimates of RIM crit- 
ical exponents. Alternatively, the non-perturbative scaling-field approach for solution of the Wilson's renormalization- 
group equation was used to study RIM critical behaviour in Ref . [Bffl . The approach similarly to e-expansion permits 
to treat model in continuous dimensions d. As a result, RIM critical exponents were found for d, 2.8 < d < 4. However, 
the scaling-field approach was not followed by more precise calculations and Ref. jSOt remains the only theoretical 
study were RIM critical exponents were obtained not on the basis of field-theoretic RG. 

Already the study of two-loop RIM massive RG functions resummed by the Chisholm-Borel procedure ||5l[] revealed 
that no difficulties connected to the degeneracy of the /3-functions are encountered. Critical exponents extracted from 
resummed 7-functions values in the fixed point were found clearly larger than the pure model ones (see table 111). 
As noted above, in the three-loop level the straightforward analysis of the /3-functions |^,^ yielded fixed point 
coordinates and critical exponents without resummation, but the obtained accuracy did not allow to estimate, for 
instance, heat capacity critical behaviour. On the other hand, the application of the Pade-Borel technique also 
encountered hardships within the three-loop level (see diagram in the Fig. [l^ ). Here, the approximant with linear 
denominator does not yield a fixed point, while another near-diagonal approximant [1/2] is unreliable because reveals 
fixed point only when processing the pole on the real axis via analytic continuation. The same expressions when 
treated by means of a Chisholm-Borel technique allowed to obtain asymptotic critical exponents (5^,^^ . 

The four-loop results were resummed both by means of Chisholm-Borel , the first confluent form of the e- 
algorithm of Wynn and Pade-Borel l55|,|6^ methods and showed close results (see tabl e ^ as well as diagram p^ . 
Whereas Pade-Borel calculations of Ref. |55f| "exploited the [3/1] Pade approximant, Ref. [|62| used a more elaborated 
generalized Pade-Borel-Leroy resummation method. The last is based on exploiting all possible Pade approximants 
in the Borel-Leroy resummation (^l])-(|6^) choosing for each of them the optimal value of the Leroy parameter q and 
then averaging the result over all values given by the approximants . 

The expressions in the five-loop approximation analyzed with the application of the flt parameter q required rather 
artificial rejecting of many approximants [|64| . For instance, using a criterion that those approximants are working, 
which provide maximal stability against variation of q, the approximant [2/ 2] was chosen to estimate u* and the 
approximant [3 / 1] to obtain v* . The analysis allowed the authors of Refs pj] to obtain the five- loop estimates for 
the RIM critical exponents (see table P^. However, the extend to the six-loop order revealed a wide gap between five- 
and six-loop fixed point coordinates |[l42(l with subsequent inconsistency of the six-loop values of critical exponents 
compared with the five-loop results of Refs [Q . This might serve as an indirect evidence of a possible non-asymptotic 
nature of the series under consideration. 

Similar resummation technique w as als o applied to the minimal subtraction RG functions (^)-(^3|) directly at 
d = 3 (scheme • 2 in the subsection VB) [^, p9| , |61 63 1. Here, asymptotic as well as effective critical behaviour was 
studied. Again, as in the massive RG case, resummation restores the presence of non-trivial fixed point in the two-loop 
approximation (see figure |l^) and preserves it in three [^ and four loops [^,^. However, in the five-loop order the 
applied Chisholm-Borel resummation scheme did not lead to a real root for the random fixed point [^ , |63| . One of 
reasons for such behaviour may be (possible) Borel-non-summability of the scries under consideration. In this case the 
four-loop approximation will be an "optimal truncation" for the resummed minimal subtraction perturbation theory 
series, similar to the non-resummed asymptotic series, e.g. in the e-expansion of 0(?i)-symmetric (p'^ model, where 
"naive" interpretation of the series truncated by e^ term leads to the best (optimal) result. 

Obtained in three- and four- [^,^,^ loop approximation values of asymptotic critical exponents are listed 
in the table [H, effective exponents were calculated [q^iIH t>y resummation of formulas (jsj), (^9|). Figure [ij shows 
the solutions u{l), v(l) of the flow equations (^5|) calculated by Chisholm-Borel resummation of four-loop /3-functions 
(|^), (H]) The flows shown in fi gure |l3| allow for several different scenarios for the values of the effective critical 
exponents (see figs. |l^, ^5|). Both in experiment as well as in computer simulations (see Tabs. ||, |l^ exponents reported 
differ and even exceed the known asymptotic theoretical values. This non-universal behaviour might be related to the 
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possible non-asymptotic behaviour found in different flows as has been suggested in Refs |3^:£j . The difference might 
be due to (i) the different temperature regions of the experiment and/or (ii) the different concentrations. The initial 
values for the couplings in the flow equations depend on the value of the concentration, especially for small dilution 
one expects the coupling v to be proportional to the concentration. If this is the case one expects a monotonous 
increase of the values of the effective exponent to the asymptotic value. A typical scenario is seen in curve 3 of Figs. 
|l|, |l4[ |l|. In this case effective exponents equal to the pure model critical exponents might be found in relatively wide 
region of temperature. Then as the attraction region of the fixed point I becomes weaker and weaker, an overshooting 
is possible and effective exponents larger their asymptotic values might be found. This scenario is predicted for larger 
dilutions and represented by the curve 6 of Figs. O, Ol U3. Curves 5, 6 correspond to situation when crossover from 



the mean field behaviour towards the random one is not influenced by the presence of a pure fixed point I. 



C. Resummation based on the conformal mapping 



The resummation procedure based on conformal mapping is widely used in the analysis of asymptotic series for 
models with one coupling, in particular for the pure d = 3 Ising model (scalar (j)'^ theory, see e.g. [126,74[ pT^ ). But 
it assumes that the behaviour of high order terms of the series expansion is known, ft is not the case for the RIM 
and hence the Pade-Borel based resummation technique was used there as described in the former subsection |Vf 
However, recently by studying the rf = dimensional quenched diluted Ising model it was shown analytically ||13q ] 
that the perturbative expansions for the free energy are Borel summable provided the summation is carried in a 
special way. The main results of Ref. ]138| | state that if the double expansion of the RIM functions in powers of u, w 
is written as 



f{u,v) = '^c^{v)u 

4=0 

oo 
3=0 



(64) 
(65) 



then the expansion (|6^) as well as expansion (|6J) at v fixed are Borel summable. These results enabled authors of 
Ref. j6^ to perform the resummation of the six-loop series (44)-(^) |]66| . Moreo ver, by no ticing that the large order 



behaviour of the series (|65|) can be derived on the basis of known asymptotics ( |60| ) [132-134] they were able to perform 
the resummation based on conformal mapping for the series in v (|65| ) to get convergent results for the coefficients 

Ci{v). 

The procedure of conformal mapping resummation of a single-variable u series is standard [126,143|. For a given 
Borel-Leroy transform S^{u) (61) the initial series may be regained from 



(66) 



Assuming the behaviour of the high order terms (|6^) one concludes that the singularity of the transformed series 
S^{u) closest to the origin is located at the point (— 1/a). Conformally mapping the u plane onto a disk of radius 1 
while leaving the origin invariant (see figure ITT), 



(1 



\l/2 



1 



(1 + au)i/2 + r 



a (1 - wy 



(67) 



substituting this into S^{u), and re-expanding in w, one receives a series defined on the disk with radius 1 in the w 
plane. This series is then re-substituted into Eq. (|66|). In order to weaken a possible singularity in t he lu-plane the 
corresponding expression is multiplied by (1 — w)^ and thus one more parameter p is introduced [ 143 |. 

However, the above procedure may be applied only to the RIM RG series for the coefficients Cj (|65|), as far as the 
large order behaviour (£^) may be derived only for series in v (|65| ) []65| . Asymptotics of the resulting series in u (^^ is 
still unknown which causes an application of the Pade-Borel resummation technique for their analysis. The last can be 
applied to the series (|65| ) as well. In summary this leads to different ways of series resummation: i) applying Pade-Borel 
method both for series (Q) and (|65|) (in Ref. this way is called a "double Pade-Borel method" ) and ii) applying 
the conformal method for series ( |65| ) and the Pade-Borel method for series ( p4| ) ("conformal Pade-Borel method"). 
Both ways were implemented in Ref. | |65| taking special attention to the choice of different Pade approximants as 
well as optimizing the results on the basis of fit parameters. In the table III we display the data for the RIM critical 
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exponents given by the authors on the basis of the analysis of the six-loop RG functions in such a way that both 
estimates obtained by the double Pade-Borel method and the conformal Pade-Borel method are included (we denote 
theme as PdBr-CM). Separately we give the results of the analysis by the double Pade-Borel method (PdBr-PdBr). 
Let us note here, that the conformal mapping technique appears to give more robust results even resumming the RG 
series for fixed v* /u* . Applying a procedure of conformal-mapping based resummation of Ref. j6^] authors of Rcf. p^ ] 
obtained for the six -lo op RG series values of critical exponents which are quite in a good agreement with the other 
estimates (see table III). However the estimates of the fixed point coordinates differ essentially in different calculation 
schemes. It leads to the conclusion |^ that probably the optimal truncation of the RIM /3-functions corresponds to 
a shorter series than a six-loop one. 



VII. CONCLUSIONS 



In the present article we have reviewed results obtained so far in the description of critical properties of a three 
dimensional weakly diluted quenched Ising model (RIM). Following the prevailing bulk of experimental, Monte-Carlo, 
and theoretical studies we focused our attention on critical exponents of the model. It can be seen from the number of 
relevant papers that the precise determination of the values of the exponents was a challenge justifying the scientific 
effort. The reason for this is two-fold, (i) The RIM allows to include in a simple model the macroscopic effects of 
disorder always present in real substances, (ii) The study of the influence of disorder on universal properties of critical 
behaviour besides practical needs is of high academic interest. It is the domain very close to critical temperature where 
even a very small amount of impurities can change the properties drastically in comparison with an ideal magnet. In 
accordance with the heuristic Harris criterion jTlj l , only when the ideal system specific heat is divergent at criticality, 
the disordered magnet is characterized by new critical exponents. The change of critical exponents of RIM is well 
established in experimental as well as Monte-Carlo and theoretical studies. However, the numerical values obtained 
show much worse self-consistency compared with the situation in studies of pure Ising magnets. This happens due 
to both technical and principal difficulties. Moreover, often the technical difficulties in determination of the RIM 
exponents are related to the principal ones! 

Though it is commonly believed that various defects are inevitable in experimental samples, the number of ex- 
perimental RIM studies is much less than those devoted to the determination of the Ising model critical exponents. 
As an explanation one might mention that the Ising model universality class includes not only magnets but also 
simple fluids, ferroelectrics, binary alloys, etc. Moreover the RIM itself can show an effective critical behaviour of the 
same type as the behaviour in the Ising model universality class (!). The data for experimentally determined RIM 
critical exponents are collected in table ||. As an peculiarity it is interesting to note that already in early studies Q| 
critical behaviour different from the pure Ising model has been observed. Surprisingly enough, since this work neither 
the experimental precision of the determination of the critical exponents has increased nor a narrower temperature 
interval around the critical point could be accessed. This can be explained by noting that starting from the end of 
eighties the researchers' attention shifted to random-field Ising model and the data for RIM were obtained only as side 
product results. It is worthwhile to note that th e the oretically predicted critical exponents for quenched disorder [ |65[ : 
V ~ 0.678,7 — 1-330 and for annealed disorder [145|: v ~ 0.708,7 — 1-391, are different but only by a small amount. 
Most probably, in real samples one encounters intermediate situations and exact coincidence between experiments on 
diluted crystalline antiferromagnets and theoretical calculations is hardly expected. 

While the existence of a new universality class for the RIM has been observed already in early experimental studies, 
this was not the case in Monte-Carlo simulations. Due to the lack of a proper finite size analysis a continuous 
dependence of the critical exponents on concentration was observed. The exponents were interpreted as effective 
critical exponents 1 68 6^ . In the following it was crucial to recognize the role of the correction-to-scaling exponent in 
the analysis of the simulation data p3|-p5[. This allowed to conclude from the numerical data on the concentration- 
independent asymptotic critical exponents . As one can note from the table |l| the improvement of computers as 
well as of the calculation algorithms allowed to increase the accuracy of critical exponents. Recently in Monte-Carlo 
simulations of the RIM attention was paid also to study of the problems of self-averaging in diluted systems [^,|8| . 

The results of experimental studies and Monte-Carlo simulations are confirmed by theoretical calculations. The 



numerical values of critical exponents obtained by different theoretical methods are collected in the table HI, Note, 



that all theoretical results were obtained within renormalization-group approach, and most of them within its field- 
theoretical formalism on the basis of the effective Landau-Ginzburg- Wilson Hamiltonian (||). Though the last approach 
encounters intrinsic obstacles such as problems with the breakdown of replica symmetry and the possible existence 
of Griffiths phase |^^, it remains the only method which can calculate the asymptotic values of exponents. For the 
pure Ising model many results are also based on other methods such as high- and low-temperature expansions. 
Another evident difference of the RIM in comparison with the Ising model lies in the applicability of different 



18 



regularisation schemes of field-theoretical renormahzation-group approach. In the case of Ising model both massive 
and minimal subtraction scheme followed by e-expansion provide consistent and reliable results. For the RIM, the 
e-expansion degenerates into -y/e-expansion 1 43 4^,^,^ which seems to be of no use for quantitative studies . 
Moreover, the initial series in the couplings, both in the massive and the minimal subtraction schemes, seem not to 
be Borel summable [135 136 13S]. Nevertheless, Padc-Borel-like resummation procedures have been employed and 
allowed to obtain convergent sequences for the critical exponents from low order series. This is evident from the 
agreement of the theoretical results obtained in different orders of perturbation theory (table |l|) with the data of 
experimental (table ^ and Monte-Carlo (table ||) studies. However, the resummation failed in higher orders which 
resulted in the conjecture that there exists an optimal truncation for the RIM RG functions [ |63[ . 

On the basis of the six-loop RG functions of the massive schemes result could only be obtained by means of a 
highly sophisticated resummation procedure |138]. In the frames of this scheme, the estimates for critical exponents are 
characterized by the same order of accuracy as for the Ising model, which are based (7|] on six-loop expansions for RG 
/3-functions and seven-loop for 7-functions. However, the determination of error-bars in the theoretical calculations 
is a difficult problem which is solved in various ways [ |62p^ . Here, the error-bars seem to measure rather uncertainty 
of specific theoretical procedures than the confidence intervals of true values. 

For a comparison of the theoretical predictions with experimental and simulation data one should use effective 
critical exponents which have been calculated for the RIM within the RG approach |^ , |63| . 

At the very end of this review we want to attract attention again to the fact that the new critical behaviour 
corresponding to the RIM universality class was experimentally observed so far only in magnetic systems like 
(anti)ferromagnets (see Table |). Still it remains a challenge to set up an experiment showing RIM-behaviour in 
other condensed matter systems. A promising example might be a fluid near its liquid-gas critical point in a porous 
medium. Recent ly, p recise experiments on the critical behaviour of liquid helium-4 near the superfluid transition in 
porous medium |146| confirmed the irrelevancy of quenched disorder. This was expected from the Harris criterion 
[147| since the specific heat exponent near the superfluid transition is roughly zero. Thus an experimental study of 
a simpl e liqu id in a porous medium may provide the first observation of a RIM critical behaviour in a non-magnetic 
system []l48|. 
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VIII. TABLES AND FIGURES 



TABLE I. The experimentally measured critical exponents of the materials, which correspond to random Ising models. The 
measurement procedures are given in the following notations: NMR — nuclear magnetic resonance; LB — linear birefringence; 
NS — neutron scattering; MS — Mossbauer spectroscopy; SMXS — synchrotron magnetic x-ray scattering; XS — x-ray 
scattering; r denotes the reduced temperature interval, where the power-law fit was done or the minimal value of the reduced 
temperature reached in an experiment. 
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TABLE IL The critical exponents of the RIM as they are obtained in MC simulations. The asterisk at concentration denotes 
that disorder was realized in a canonical manner. 
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TABLE in. The theoretical values for asymptotic critical exponents of the RIM. nLA denotes the nth order in loopwise 
approximation within massive ('mass') and 3d minimal subtraction ('MS') schemes of field-theoretical renormalization group 
approach. The resummation procedures are given in the following notations: ChBr - Chisholm-Borel; PdBr - Pade-Borel; 
AW - e algorithm of Wynn, CM - Borel transformation with conformal mapping. SF stands for Golner-Riedel scaling field 
approach, superscript at the correction-to-scaling exponent uj indicates that the real part of corresponding complex number 
is shown. 
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FIG. 1. Weakly diluted quenched 3d Ising model (random Ising model: RIM) describing a system of scalar "spins" randomly 
distributed in sites of a three-dimensional cubic lattice and fixed in certain positions. 




FIG. 2. RIM fixed points qualitative structure. Gaussian fixed point G is stable for d > 4, stable fixed point P can not be 
reached from the initial coupling values u > 0, v < (this region as well as other unphysical for RIM regions are shown in grey 
tone), pure Ising fixed point I is unstable whereas random fixed point R is both stable and accessible (stable fixed points are 
shown by boxes). 
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FIG. 3. Neutron scattering measurements of the susceptibility x(0) in Mno.75Zno.25F2 . The sohd hues are the resuhs of 
fits to single power laws with exponents 7 ~ 1.364 above and below Neel temperature Tjv. The critical behaviour is governed 



by RIM asymptotic critical exponents over the reduced temperature interval 4-10 < [t[ < 2 ■ 10 
the Ref. i. 



The figure is taken from 
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FIG. 4. Neutron scattering measurements of the longitudinal inverse correlation length k in Mno.5Zno.5F2 above and below 
Tn The solid lines are the results of fits to laws with the exponents v' ~ 0.76 for T < Tn and ~ 0.75 for T > Tjv. The 
data give v' = v within the errors (see Table ^. The law holds in the reduced temperature interval 5 • 10~^ < |r| < 10~^. The 
asymptotic region still is not reached as the observed effective exponents values show. The figure is taken from the Ref. [^. 
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FIG. 5. The temperature derivative of the optical birefringence of Mno.55Zno.45F2 at zero magnetic field. In the main 
part of the figure d(An)/dT is shown for the critical region (indicated by box in the inset). The measurements bring about a 
cusp- like behaviour of the RIM specific heat with the critical exponent a — —0.09 ± 0.03. The figure is taken from Ref. ]13|. 
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FIG. 6. Obtained in MC simulations |29| RIM effective critical exponents /3 (dots) and 7 (squares) for different dilutions 
p. All exponents change continuously with dilution and are clearly different from their pure system values. The figure is taken 
from Ref. H. 
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FIG. 7. Determination of the correction-to-scaling exponent uj in the MC simulations of Ref. |33|. Quantity Qa^j — 2^+^/" 
is plotted for different dilutions 0.4 < p < 1 and lattice sizes 8 < L < 128. One can see that at p = 0.9 the system is crossing 
over from the pure Ising fixed point to the diluted one even for L — 128. The solid lines correspond to a fit a; = 0.37 yielding 
the same infinite volume extrapolation for all p < 0.8. The figure is taken from Ref. [B3|. 
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FIG. 8. Two-loop integrals of the RG functions in massive scheme (formulas (fl8[)-(|50[)) as functions of space dimension d, 
Ref. M. 
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FIG. 10. The lines of zeros of non-resummed (left-hand column) and resummed by the Chisholm-Borel method (right-hand 
column) RIM massive /3-functions in different orders of the perturbation theory: one- and two-loop approximations. Circles 
correspond to Pu = 0, thick lines depict /3„ = 0. Thin solid and dashed lines show zeros of the analytically continued functions 
l3u and Pv respectively. One can see the appearance of the mixed fixed point m > 0, ?; < in the two-loop approximation for 
the resummed /3-functions. The figure is taken from Ref. [pO[. 
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FIG. 11. The lines of zeros of non-resummed (left-hand column) and resummed by the Chisholm-Borel method (right-hand 
column) RIM massive /^-functions in three- and four-loop approximations. The notations are the same as in the figure hd. Close 
to the mixed fixed point the behaviour of the resummed functions remains alike with the increase of the order of approximation. 
This is not the case for non-resummed functions. The figure is taken from Ref. ||6c| . 
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FIG. 12. Pade-Borel resummation technique applied to the /3-functions of the figure ^. Resummation restores the presence 
of the fixe d po int u* ^ 0, v* = and results in the appearance of a new stable fixed point u* 7^ 0, v* ^ 0. The figure is taken 
from Ref. [129]. Note that the authors of Ref. [12£| exploited different in comparison to formulas 
couplings 3/2u and v —4 v. 
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FIG. 13. Flow lines for the RIM. Fixed points G, I are unstable, fixed point R is the stable one. The figure is taken from 
Ref. [^. Note that the authors of Ref. exploited different in comparison to formulas (^o|)-(^) normalization of couplings 
11 — > 3 u and « — > 8/3 ii. 
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FIG. 14. Effective exponent v^jj versus logarithm of the flow parameter I for the flows shown in Fig. |l3|. The figure is 
taken from Ref. [ psf . 
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FIG. 15. Effective exponent 7e// versus logarithm of the flow parameter £ for the flows shown in Fig 
taken from Ref. Ipsf . 
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FIG. 16. The fecundity of Pade-Borel resummation for the RIM /3-functions in massive scheme. The numbers of a row M 
and of a column N correspond to a Fade approximant [M/N] representing (M+N)-loop approximation. Symbol "-" depicts 
that no fixed point values are available, "±" and "~" show that a method either uses analytic continuation due to poles on the 
real semi-axis or yields values very far from expected. The fruitful approximants are shown by dark colour and/or appropriate 
citations. Though working approximants are distributed rather stochastically, numerical results on their basis show evident 
convergence (table ffl|). 
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FIG. 17. Conformal mapping of the cut-plane onto a disc leaving the origin invariant, as prescribed by formulas (|67|). See 
the text for details. The figure is taken from Ref. [144]. 



37 



